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Abstract. Gravitational waves are radiative solutions of space-time dynamics predicted by Einstein's theory of General 
Relativity. A world-wide array of large-scale and highly-sensitive interferometric detectors constantly scrutinizes the geometry 
of the local space-time with the hope to detect deviations that would signal an impinging gravitational wave from a remote 
astrophysical source. Finding the rare and weak signature of gravitational waves buried in non-stationary and non-Gaussian 
instrument noise is a particularly challenging problem. We will give an overview of the data-analysis techniques and associated 
observational results obtained so far by Virgo (in Europe) and LIGO (in the US), along with the prospects offered by the up- 
coming advanced versions of those detectors. 
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Einstein's theory of General Relativity introduces the 
concept of a deformable and evolving space-time. The 
dynamics of space-time is prescribed by the Einstein 
equation. In the linearized gravity which assumes small 
deformations in a nearly flat space-time, this equation re- 
duces to the wave equation which therefore evidences the 
existence of radiative solutions. The latter are referred 
to as gravitational waves (GW) and can be phenomeno- 
logically seen as propagating disturbances of space-time 
itself. The theory also predicts that GW are transverse 
waves, that they nominally propagate at the speed of light 
and possess two independent polarizations [1,2]. 

GW have never been directly detected i.e., through the 
measurement of their effect on a man-made instrument. 
Strong evidence of their existence has been provided by 
the observation of the famous Hulse-Taylor pulsar binary 
(PSR B1913-H16) [3]. The decay rate of the binary orbital 
period is in remarkable agreement with the predicted 
evolution obtained under the assumption that this system 
radiates energy away in the form of GW. 

The direct search for GW made notable progresses 
with the advent of dedicated instruments based on high- 
precision laser interferometry such as LIGO and Virgo 
(see [4, 5] for a detailed review). With the ongoing instal- 
lation of a new and improved generation of those instru- 
ments, the first discovery is expected within the decade. 

While electromagnetic waves are produced by accel- 
erated charges, GW are produced by accelerated masses. 
Very large masses and relativistic velocities are neces- 
sary to generate GW at a detectable level. For this reason, 
the current projects aiming at detecting GW target poten- 
tial astrophysical sources involving very dense and com- 
pact objects such as neutron stars or black holes. Very 



energetic astrophysical events such as the coalescence of 
neutron star and/or black holes binaries, or stellar core 
collapses are expected to be the source of intense and 
short-duration bursts of GW [4]. 

Because of the limited rate of occurrence for such 
events, searching for such transient GW in the LIGO and 
Virgo data essentially consists in searching for rare and 
weak signals at the detectability limit. This article gives 
the state-of-the-art of the search for GW transient sig- 
nals with a focus on the related data analysis challenges. 
Searches for long-lived signals such as periodic GWs 
from rotating neutron stars, and stochastic GW back- 
grounds are beyond the scope of this paper. We first give 
some introductory material with a general presentation of 
the detectors in Sec. 1 and a review of the relevant astro- 
physical sources in Sec. 2. Sec. 3 gives an overview the 
major problems faced when searching for transient GW 
along with the data analysis methods deployed to address 
them. 



1. INTERFEROMETRIC GW 
DETECTORS 

The first generation of interferometric GW detec- 
tors comprises five large-scale instruments in to- 
tal (see Fig. 1). The US-based Laser Interferometer 
Gravitational- Wave Observatory (LIGO) [6] includes 
three kilometric-scale instruments located in Livingston, 
Louisiana (labelled LI) and Hanford, Washington (the 
latter hosting two interferometers in the same vacuum 
enclosure with labels HI and H2). The French-Italian 
project Virgo [7] has one instrument of the same class 




FIGURE 1. Geographic location of the current and future 
GW interferometric detectors. This world map displays the 
location of the four sites of the first generation detectors (LIGO 
H and L, Virgo and GEO), and six sites of the second generation 
(complemented by LIGO I and KAGRA). The future detector 
LIGO India is still pending approval and its exact location is 
yet to be determined. Credits: [9] 

located in Cascina near Pisa, Italy (labelled VI). This 
set of kilometer- scale instruments is complemented by a 
detector with more modest dimensions (several hundreds 
of meters): GEO [8] (labelled Gl), a German-British 
detector in operation near Hanover, Germany. 

Despite major differences in the technologies in 
use, all those instruments measure gravitational waves 
through the same principle. They all sense the strain that 
a passing GW exerts on space-time by monitoring the 
differential length 5£ of the optical path followed by 
two laser beams propagating along orthogonal directions 
over a distance L. This is performed by letting the two 
beams used for the measurement interfere similarly to 
the Michelson-Morley experiment. The interference is 
closely related to the difference in the phase accumulated 
by the two beams before they combine and hence to 
the difference in their optical paths. The measurement 
of the interference light power allows that of 5i with 
high accuracy. Measurement noises (mainly the ther- 
mal noise due to the Brownian agitation of the atoms 
constitutive of the optics and the shot noise due to the 
quantum nature of light) can be reduced to reach the 
level ofh = 51 /L ~ 10~^^ where h is detector response, 
homogeneous to the a-dimensional GW amplitude. The 
best sensitivity is achieved in a frequency band ranging 
from ~ 100 Hz to 1 kHz approximately (see Eig. 2 - 
bottom). 

The detector response is a linear mixture h = F^h^ -\- 
Fxhx of the two GW polarizations /z+ and hx- The an- 
tenna pattern factors and F^ characterize the way the 
wave polarizations couple to the detector. The coupling 
^ — {F^ -\- F^y^^ < 1 is maximum for waves coming 
from the normal to the detector plane and is minimum 
(and exactly zero) for waves from the four "blind" di- 
rections associated to the two bisectors of the detector 
arms. GW detectors are non-directional instruments as 
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FIGURE 2. (top) Time line of the data takings completed 
so far. Credits: [10] (bottom) Sensitivity achieved by LIGO 
and Virgo during their last science data taking S6A^SR2-3 

[11]. 

^ > 1 /2 for more than half of the sky. 

The first generation detectors have conducted a series 
of science data takings reaching an integrated observa- 
tion time of about 2 years (see Fig. 2 - top). The data 
takings are coordinated in order to maximize the obser- 
vation time with the three most sensitive detectors oper- 
ating while always maintaining at least one detector in 
"astro- watch" mode in case of an outstanding galactic 
event. 

The first generation of detectors has now been decom- 
missioned and it is currently being replaced by a second 
generation of "advanced" detectors. Thanks to major up- 
grades in their infrastructure and instrumentation, a ten- 
fold increase in sensitivity is expected with the advanced 
detectors as indicated in Fig. 3. The GW amplitude de- 
caying inversely with the distance, this corresponds to a 
factor of thousand in the observable volume of Universe 
and hence in the number of detectable sources. Advanced 
detectors are likely to detect several tens and possible 
several hundreds of sources as we will see in the next 
Section. 

The installation of the advanced LIGO detectors [12] 
should be completed by the end of 2013 and a first sci- 
ence run is likely to take place in 2015. The original plan 
was to install two four-kilometer detectors at Hanford 
site, but there is now a proposal (still to be approved by 
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FIGURE 3. Projected sensitivities for the advanced LIGO 
and advanced Virgo detectors compared with the design 
sensitivity of their initial version. 



U.S. and Indian institutions) to move one of those de- 
tectors at a new observatory in India. If this plan mate- 
riaUzes, the third detector at the Indian site would start 
operation around 2020. Advanced Virgo plans to have 
a robustly operating detector in 2015 and to begin col- 
lecting science data as soon as possible after that [13]. 
GEO foresees a program of upgrades "GEO-HF" [14] 
which focuses on improving the detector sensitivity at 
high-frequencies thanks to a larger laser power and the 
use of "squeezed light". The network of advanced detec- 
tors will be completed by the Japanese KAGRA detec- 
tor [15] which has the specificity of being installed un- 
derground in the Kamioka mine (where the seismic mo- 
tion is much lower than at the surface) and to operate at 
cryogenic temperature to reduce thermal noise. An ini- 
tial three-kilometer room-temperature interferometer is 
expected to be operational by 2015, with the full cryo- 
genic interferometer ready to start taking data by 2018. 

2. ASTROPHYSICAL SOURCES OF GW 
TRANSIENTS 

Phenomenologically, GW emission arises from relativis- 
tic bulk motion. At lowest order, GW can be related to 
variations in the quadrupolar moment of the mass distri- 
bution [1]. Therefore, GW sources have to present some 
degree of non-axisymmetry. In this section, we review 
the astrophysical scenarios giving rise to GW emission. 
The coalescence of neutron-star and/or black-hole bina- 
ries similar to the Hulse-Taylor binary mentioned previ- 
ously is often considered the most promising one. 

The last minutes before the system merges give rise to 
the emission of an intense burst of GW. Post-Newtonian 



expansions of the binary dynamics [16, 17] are used to 
predict the gravitational waveforms radiated during the 
inspiral phase which precedes the merger. The GW sig- 
nature consists in a chirp signal whose frequency sweeps 
toward high values according to a power law at first order 
. A substantial amount of energy is radiated in the follow- 
ing phase when the two bodies merges into a black-hole. 
In this highly relativistic phase, the perturbative treat- 
ment of binary dynamics is not valid anymore and one 
has to resort to numerical simulations. The process is 
concluded by the ring-down phase during which the re- 
sulting distorted black-hole radiates away its asymmetry 
down to equilibrium. During the whole coalescence pro- 
cess, an stellar-mass binary with equal masses radiates 
away of order of a percent of its rest mass [4]. 

Although binary systems are fairly common, only a 
small fraction eventually forms a compact binary that is 
sufficiently tight to coalesce in less than Hubble time. 
A survey of population estimates [18] gives a "realis- 
tic" rate of one neutron- star-neutron- star coalescence^ 
per 10,000 year per galaxy equivalent in size to the Milky 
Way. GW detectors can ideally observed those binary 
systems up to a distance of ~ 30 Mpc and ~ 440 Mpc for 
initial and advanced detectors resp. [18]. Converted into 
a rate of detectable coalescences, this leads to ~ 0.02 
events per year with the first generation of (initial) de- 
tectors and to ~ 40 events for the second generation (ad- 
vanced). Large error bars are attached to those estimates 
reflecting the weakness of the observation constraints we 
have about those systems. The above stated rates can then 
be 10 times smaller or larger in the "pessimistic" or "op- 
timistic" scenarios respectively. The "realistic" rates pre- 
sented above are corroborated by the ones derived as- 
suming that compact binary mergers are the progenitors 
of short-hard gamma-ray bursts (GRB) [18]. 

Gravitational stellar-core collapse is another potential 
source of GW if some degree of non-axisymmetry is ex- 
hibited during this process. The simulations required to 
make reliable predictions of the emission levels are very 
challenging as they have to incorporate many physical in- 
gredients including relativistic magneto-hydrodynamics 
and a detailed treatment of neutrino transport and nu- 
clear iterations [19]. The current realistic estimate of the 
amount of radiated GW energy is of order lO^^M© and 
corresponds approximately [20] to a distance reach of or- 
der ~ 10 kpc with the initial detectors, ~ 100 kpc with 
the advanced detectors. The detectable sources are there- 
fore located in the Galaxy. 

Another potential source of GW bursts are "neutron- 
star quakes" [21] during which the vibrational normal 
modes of a neutron- star are excited and damped by GW 



Similar rates are obtained for the other types of systems mixing 
neutron-stars and/or black-holes. 



emission. Star quakes may origin from the disruption of 
the star crust due to the sudden rearrangement of the 
magnetic field of a highly-magnetized neutron- star (mag- 
netar). Cosmic string cusps may be also listed among the 
potential GW burst sources [22]. 

3. SEARCHES FOR GW TRANSIENT 
SIGNALS 

3.1. Time-series analysis 

In detection problems, the availability of a priori in- 
formation plays a major role. We have seen in Sec. 2 
that the GW signature from coalescing binaries of neu- 
tron stars and/or black holes have a specific time evo- 
lution which can be predicted with good accuracy. This 
morphological information helps to distinguish a real 
GW signal from the instrumental or environmental noise. 
The search for known signals is efficiently performed by 
matched filtering techniques [23] which cross-correlates 
the data with the expected "template" waveforms ob- 
tained from the source model. 

Because of the highly-relativistic dynamics associated 
with the production of GW, some of the expected GW 
waveforms are difficult to predict with accuracy. This 
calls for detection methods that are robust to the model 
uncertainties. Excess power methods essentially consist 
in searching for a broad family of GW waveforms by 
scanning a time-frequency map for transient excursions. 
The time-frequency map is obtained by projecting the 
data onto a dictionary of elementary waveforms that 
tiles the time-frequency plane. Several types of dictio- 
nary have been tested including local cosines [24], sine- 
Gaussian wavelets [25], orthogonal wavelet packet bases 
[26] or chirplets [27]. Real GW signals are unlikely to 
correlate exactly with one element in the dictionary, but 
with several of them. Clustering algorithms are generally 
applied to harvest the signal energy scattered over several 
elements [28, 29, 30]. 

The time-frequency dictionaries mentioned above are 
composed of "generic" elementary waveforms mainly 
motivated by mathematical or algorithmic arguments. 
Astrophysically motivated dictionaries can be obtained 
by extracting the relevant information from catalogs of 
GW signals developed through numerical simulations 
[31,32, 33, 34]. 



3.2. IMulti-detector analysis 

We described the basic ideas employed to analyze the 
data stream from individual detector. A gain in sensitivity 
is expected from the availability of joint observation by 



multiple detectors. This section discusses several aspects 
related to the combined analysis of multiple detector 
data. 



3.2.1. Coherent analysis of multiple data streams 

We already mentioned that the detector receives a mix- 
ture of both GW polarizations which depends on the rel- 
ative orientation and alignment of the detector and wave. 
Since the detectors are not co-planar and co-aligned, they 
couple differently to the incoming wave resulting in ob- 
served responses with different initial phases and am- 
plitudes. Because of its finite speed, GW reaches the 
detectors at different times. All those differences can 
be exploited using coherent analysis techniques to im- 
prove the overall sensitivity. Those techniques consist 
in compensating the phase shift and delay of the vari- 
ous responses to align them in time and phase assum- 
ing a given direction-of- arrival. The resulting data stream 
are combined so that the sum operates constructively for 
GW signals from the selected direction. The data stream 
which results from the coherent combination maximizes 
the signal-to-noise ratio (SNR). The combined stream 
can be then analyzed using methods inspired from the 
single detector case, i.e., excess power methods for the 
unmodelled GW bursts [35, 36] and matched filtering 
techniques [37, 38] for inspiralling binaries. The coher- 
ent analysis being directional (each coherently combined 
stream is associated to a given direction), the outcome is 
a probability (pseudo-)distribution over the sky usually 
referred to sky map from which the most likely location 
of the source can be extracted. 



3.2.2. Mitigation of non-Gaussian/non-stationary noise 

The noise of the real instruments is far from the ideal 
properties of stationarity and Gaussianity we expect from 
the main fundamental (thermal and quantum) noises. The 
tails of the noise distribution is dominated by a non- 
Gaussian and non- stationary component consisting in a 
large number of transient noise excursions commonly 
called glitches. Glitches are produced by a variety of en- 
vironmental and instrumental processes, such as upcon- 
version of seismic noise or saturations in feedback con- 
trol systems. Since glitches occasionally occur nearly si- 
multaneously in separate detectors by chance, they can 
mimic a gravitational- wave signal. 

The population of glitches is difficult to model. The 
size and the large complexity of GW detectors makes 
this modelling even more difficult. GW detectors being 
extended instruments over kilometers, it is hard to com- 
pletely isolate them from the outside world and the sur- 
rounding anthropic activity. Therefore, the accurate mod- 



elling of the non-Gaussian/non-stationary noise back- 
ground is for now out-of-reach. It has to be mitigated and 
this can be done at least partially by using multiple data. 

It is possible to calculate combinations of the data 
from multiple detectors where the GW signals from all 
detectors interfere destructively in the sum. The GW 
signal thus cancels, but not background glitches. The 
energy in these "null" stream(s) may be used to reject 
or down- weight events not consistent with a gravitational 
wave [39, 36]. The success of such tests depend critically 
on having several independent detectors of comparable 
sensitivity. 



3.2.3. Use of multiple detector data for background 
estimation 

We explained earlier that the accurate modelling of 
non-Gaussian non- stationary noise is out-of-reach. The 
remaining part of the glitches that cannot be identified by 
the coherent techniques described in the previous section 
constitutes the dominating background noise in burst 
searches. This background has to be estimated. However, 
GW signals cannot be turned off: the detectors cannot 
be shielded from them. Therefore, we don't have "noise- 
only" data at disposal for background estimation. 

Nevertheless, the background can be estimated thanks 
to the availability of multiple data streams by time shift- 
ing one detector's data. The time- shift is chosen to be 
much longer than the time-of- flight between detectors (~ 
30 ms) and coherence time scale of the detector noise 
seconds). Time-shifted (or "time slide") analysis leave 
only triggers due to accidental coincidences of instru- 
mental glitches. The contribution from real GW signals is 
practically erased. By repeating the analysis many times 
with different time-shifts, we get an accurate estimate of 
the rate of background events. For sufficiently large time 
shifts, each trial can be considered independent to the 
other. However, the number of time- slides cannot be in- 
creased indefinitely as a significant correlation between 
time-slides will occur above a certain level [40]. 

The p- value measuring the significance of a GW event 
can be computed by computing the fraction of louder 
background events from the time-slide analysis. 
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FIGURE 4. Examples of background distributions for the 
coherent Wave-Burst algorithm [36] for a three-detector 
network (LIGO-H and L and Virgo) during S6-VSR2/3 
data taking. This distribution is given as a function of the cor- 
related amplitude p homogeneous to the signal-to-noise ratio. 
Hatch area: background before any glitch rejection scheme 
is applied. Black area: after the "null-stream" glitch rejec- 
tion (see Sec. 3.2.2). Gray area: after data-quality flags (see 
Sec. 3.3). CAT2 and 3 refers to the different categories of the 
data-quality flags whose description goes beyond the scope of 
this article. Bold curve: expectation if noise is stationary and 
Gaussian. Credits: [41]. 



GW channel with the auxiliary channels can help de- 
termine the origin of noise artifacts and how the orig- 
inal disturbance couples into the detector [41]. A sig- 
nificant number of noise sources are identified a poste- 
riori after the science data taking is done. Those noise 
sources cannot be mitigated by fixing the instrument. 
Instead, this leads to the development of a data-quality 
flag which when "raised" indicates that the data are im- 
proper, and any event occurring at that time should be ve- 
toed. This provides also an important resource for back- 
ground glitch rejection. Data-quality flags with a large 
1) efficiency (percentage of glitches vetoed by the 
flag) over dead-time (fraction of science time rejected by 
the flag) ratio of are particular interest [41]. About 200 
data-quality flags are used in GW burst searches. Fig- 
ure 4 shows the background improvement after vetoing. 



3.3. Data quality 

Besides the gravitational- wave channel h{t), hundreds 
of auxiliary channels including microphones, seismome- 
ters, magnetometers, etc. are recorded at any given time 
during science data takings. Those channels can be used 
to get an image of the operational and environmental sta- 
tus of the detector. The observed correlation between the 



3.4. Results 

LIGO and Virgo conducted two joint science runs, la- 
belled 85 for LIGO and VSRl for Virgo for the first run, 
S6 and VSR2/3 for the second. A total of T = 635 days of 
observing time have been analyzed [43, 44]. No GW de- 
tection has been claimed yet. Upper limits (at 50% con- 
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FIGURE 5. Cumulative event rate associated to the 
search for compact binary coalescences. The "candidate" at 
pc ~ 12.5 is a simulated signal inserted as part of a "blind 
injection challenge exercise" (see text). Credits: [42] 



fidence level) on the GW strain obtained from an all- sky 
all-time GW burst search have been set. It is slightly be- 
low hrss < 5 X 10~^^ Hz~^/^ for frequency at about 200 
Hz, where the bound is on the root-mean- square ampli- 
tude h^^^ jdth\(t) ^h\(t) of the two GW polariza- 
tions at Earth. While the exact result depends on the as- 
sumed GW model (here, a generic sine Gaussian wave- 
form characterized by its central frequency), it remains 
comparable for waveforms with a similar spectrum. 

This upper limit placed on a local quantity {hyss at 
Earth) can be translated into astrophysical constraints. 
For instance, lower limits on the radiated energy Eqw 
by generic sources of linearly polarized GW located at 
distance d. Averaging over the source inclination, the 
above strain limit corresponds to Eqw = 2 x IO'^MqC^ 
for galactic sources at distance J = 10 kpc, and Eqw = 
5 X 1O~^M0C^ for source in the Virgo cluster with J = 15 
Mpc. Those estimates are comparable to the expected 
GW-radiated energy from core collapses and mergers of 
stellar-mass compact objects respectively. 

The same data, when searched specifically for inspi- 
ralling binaries of neutron stars, leads to an upper limit 
on the rate of such astrophysical events of ^90% = 1.3 x 
10~^yr~^Mpc~^ [45, 42] which is still two orders of 
magnitude larger than the rate estimate obtained from 
population models [18]. Fig. 5 shows the cumulative 
rate of events detected by the matched-filtering proce- 
dure outlined in Sec. 3.1 in coincidence in the HI and LI 
detectors during four months of S6A^SR2-3 data taking. 
This distribution is displayed as a function of the rank- 
ing statistic pc which combines the signal-to-noise ratios 
measured at both detectors. The distribution of candidate 
events (triangle) is superimposed to a background esti- 
mate (black dots) with error bars (in gray). The last tri- 
angle on the right-hand side of the plot at pc ~ 12.5 is 
a candidate event detected with false-alarm rate of 1 in 
7,000 years [42]. It was known in advance that a small 
number of fake GW signals might be added "blindly" to 



the data. The exact time and characteristics of these sig- 
nals were only known by a small group of people sworn 
to silence until the eventual "opening of the envelope". 
The envelope was not empty and contained the detected 
event code-named GW 1009 16 [46]. Thanks to this ex- 
ercise, it was possible to test the entire decision-making 
chain up through the preparation of a publication for re- 
lease. 



4. ONLINE ANALYSIS AND RAPID 
ELECTROMAGNETIC FOLLOW-UP 

Sources of GW are likely sources of other kinds of emis- 
sions, such as electromagnetic waves or jets of high- 
energy particles. The possible connection between com- 
pact binary coalescences and GRB is an example [47]. 
This motivates cross-correlating GW with other types of 
observations in the electromagnetic or neutrino spectra 
see e.g., [48, 49] for recent results. We will briefly report 
here on rapid follow-up observations seeking electro- 
magnetic counterparts to GW candidate events. A low- 
latency analysis pipeline was operated for the first time 
during the last data taking [50]. It allows to generate 
alerts "on the fly" within 20 minutes of the associated 
GW candidate event. IVLajor changes were required with 
respect to the original off-line pipelines. The most prob- 
able direction of the source along with an error box were 
communicated to a dozen of partner observatories [50] 
including radio-telescopes, wide-field optical telescopes 
and the X/gamma-rays satellite Swift. This has led to 
follow-up observations which have been scanned for 
transient excursions. This exercise has been extremely 
useful and will help to prepare the exciting future of 
multi-messenger astronomy with the advanced detectors 
[51,9]. 
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